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SPECIMENS BY BOUNDARY COLLOCATION 
by Bernard Gross and John E. Srawley 
Lewis Research Center 

SUMMARY 

A boundary- value- collocation procedure was applied to the Williams stress function 
to determine values of the stress-intensity factor K for single-edge cracks in 
rectangular- section specimens subjected to three-point bending. The results are pre- 
sented in terms of the dimensionless quantity = K^B^W^/M^ where B and W are 
the specimen thickness and depth and M is the bending moment at midspan. The values 
of Y Z as a function of relative crack depth a/W for three-point bending are appre- 
ciably lower than the corresponding values for pure bending (determined previously by 
the same method) and decrease as the ratio of support span to specimen depth S/W 
decreases. Plots of Y against a/W are given for values of a/W up to 0. 5 and S/W 
equal to 4 and 8. 

The results were relatively insensitive to variations in the spread of the midspan 

load contact region, which was assumed to be related to the yield strength of the 

material. The results agreed fairly well with published results derived from experi- 

2 

mental compliance measurements; one set gave higher values of Y than the present 
method, and the other set gave lower values. The plane- strain fracture toughness meas- 
urement capacity of three-point bend specimens is somewhat lower than that of four- 
point bend specimens, but the difference is of negligible practical importance. 


INTRODUCTION 

Various types of crack-notch specimens are used for Kj c plane-strain fracture 
toughness testing of materials (ref. 1). (Symbols are defined in appendix A. ) The value 
of Kj c determined in such a test is equal to that value of the crack-tip stress-intensity 
factor K at which the crack becomes unstable and extends abruptly in the opening mode. 
The stress- intensity factor K is proportional to the applied load and is a function of the 
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(a) Four-point bending (regarded as practically equivalent 
to pure bending in central test region near crack). 
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(b) Three-point bending (conventional force diagram). 
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specimen dimensions, particularly those of the 
crack, which depends upon the configuration of the 
specimen and manner of loading. The determi- 
nation of an expression for K in terms of these 
factors for a particular type of specimen will be 
referred to as the K calibration for that speci- 
men type. Various experimental and analytical 
methods of stress analysis for K calibration 
have been developed (ref. 2). 

One important class of Kj c specimens, 
which may be referred to briefly as single- edge- 
notch specimens, includes those of rectangular 
cross section having a single crack-notch extend- 
ing from one edge. The different types of speci- 
men within this class are distinguished primarily 
by the manner of loading. Previous reports by 
Gross, Srawley, and Brown dealt with the appli- 
cation of the boundary collocation method of 
stress analysis for K calibration to single- edge- 
notch specimens loaded in uniform tension 
(ref. 3) and to such specimens loaded either in 
pure bending or in combined bending and tension (ref. 4). 

In practice, single- edge- notch specimens are often tested in three-point bending, 
whereas the results for pure bending given in reference 4 apply only to ideal four-point 
bending. In four- point bending (fig. 1(a)), providing that the loads are applied at positions 
sufficiently distant from the crack, the stress-intensity factor K depends upon four 
variables only, namely, the applied bending moment M, the crack length a, the specimen 
depth W, and the thickness B. The K calibration is conveniently expressed in terms of 
the dimensionless quantity = K^B^W^/M^, which is a function of a/W only. 

In three-point bending (fig. 1(b)), however, there are two additional independent 

variables which might affect K. One of these variables is the support span S. For a 

given value of the maximum bending moment M = PS/4, there is a bending moment 

gradient and a shearing force, both inversely proportional to S. Thus, in terms of the 

o 

dimensionless form of the K calibration, it is to be expected that Y will be a function 
of S/W as well as a/W. The second additional variable concerns the distribution of 
contact pressure around the nominal position of the central loading point. The load is 
usually applied through a hard, cylindrical roller, and for the present purpose it is 
assumed that the specimen behaves in the contact region as a rigid-plastic solid with a 
well-defined yield strength Oyg. The applied load P can then be regarded as evenly 


(c) Three-point bending (distributions of boundary stresses 
in central region extending distance V on either side of 
crack). 

Figure 1. - Diagrammatic representation of beam loading. 
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TABLE I. - COMPARISON OF RESULTS OF PRESENT WORK WITH 


CORRESPONDING RESULTS FROM REFERENCES 4, 5, 6, AND 7 
IN TERMS OF DIMENSIONLESS QUANTITY E^bV/M 2 AS 
FUNCTION OF RELATIVE CRACK LENGTH 


Relative 

Results of 

crack 






length, 

Boundary 

Reference 4 

Reference 5 

Reference 6 

Reference 7 

a/W 

collocation 






(a) 

(b) 

(a) 

(c) 




Y 

2 = k^Wm 2 


0.10 

11.70 

12.40 

9.44 

10.88 

10. 08 

.15 

17.30 

18. 50 

13.92 

17.28 


.20 

23.47 

25.30 

19.20 

24. 48 

22.29 

.25 

30.78 

33.20 

26.08 

33.12 


.30 

39.84 

42.80 

35. 20 

41.92 


.35 

51. 54 

55. 20 

48.16 

57.28 


.40 

67.18 

71.41 

63.04 



.45 

88.78 

92.70 

83.36 



. 50 

119.74 

123.01 

108.80 


140.48 


a Ratio of support span to specimen depth, 8. 
kpure bending. 

c Ratio of support span to specimen depth, 10. 


distributed over a rectangular area of breadth B (the thickness of the beam) and length 
2y o (fig. 1(c)), such that 2y Q B is equal to P/cr y g. Hence, the third variable which af- 
fects the value of can be taken as 2y Q /W. 

Calibrations for K were computed for S/W values of 4, 6, 8, and 10 to cover 
what was considered to be a practical range. In the past, a ratio of 8 has been commonly 
employed, but there is an increasing tendency toward the use of smaller ratios in order 
to conserve test material. On the other hand, the accuracy of K calibrations for S/W 
less than 4 was considered dubious because of the increasing difficulty of representing 
the physical loading conditions accurately as this ratio becomes smaller. 

The ratio 2y Q /W has an upper limit which is related to the condition for a valid 
K lc measurement, namely, that the nominal stress at the crack tip should not exceed 
the yield strength of the material (ref. 1). From this condition it can be shown that the 
upper limit of 2y Q /W is equal to 2(1 - a/W)^W/3 S. Calibrations for K were conducted 
in parallel by using this upper limit and taking zero to be the lower limit. The differences 
between the results for the two limits were so small that it was not necessary to consider 
intermediate values of 2y Q /W. 
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Results of experimental compliance measurements on three- point bend specimens 
have been published by Irwin, Kies, and Smith (ref. 5) and by Kies, Smith, Romine, and 
Bernstein (ref. 6). Limited results of an analytical study by H. F. Bueckner have been 
published by Wundt (ref. 7). As shown in table I, there is sufficient lack of agreement 
between K calibration values derived from the results of these three references to war- 
rant the undertaking of the present study. Furthermore, these references provide no 
information about the extent to which the K calibration depends on the parameters S/W 
and 2y 0 /W. In the interest of accurate measurement of Kj c , it is important that the 
extent of the influence of these parameters should be known. 

ANALYTICAL AND COMPUTATIONAL PROCEDURE 


The method of analysis consists in finding a stress function x that satisfies the bi- 
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harmonic equation V x = 0 and also the boundary conditions at a finite number of sta- 
tions along the boundary of a single- edge-notched specimen, such as shown in figure 1 
(p. 2). The biharmonic equation and the boundary conditions along the crack are satis- 
fied by the Williams stress function (ref. 8). Because of symmetry (fig. 1), the coef- 
ficient of the sine terms in the general stress function must be zero; hence, 
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The stresses in terms of x obtained by partial differentiation are as follows: 
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The boundary collocation procedure consists in solving a set of 2m simultaneous 
algebraic equations which correspond to the known values of x and either 9x/9x or 
9x/3y at m selected stations along the boundary ABCC’D of figure 1(c); thus, values for 
the first 2m coefficients of the Williams stress function are obtained when the remaining 
terms are neglected. Only the value of the first coefficient d^ is needed for the present 
purpose since the stress-intensity factor K is equal to - yihrdp as shown in reference 4. 

The required values of x and its first derivatives at the m selected boundary 
stations were obtained from distributions of bending moment, shear, and contact 
stresses (fig. 1(c)), equivalent to the concentrated loads of the conventional force diagram 
(fig. 1(b)). The equations for these boundary values, in dimensionless form, areas 
follows: 


Along AB 


Along BC 


X = 0; W 8 -{x/?). = 0 

P 9x 


^ - v) — (W - 2a) + — (W ~ a) a\1 
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The distance V of the end boundary from the crack (fig. 1(c)), was chosen to be 
approximately 1. 5 W, the exact value being different for different values of a/W as a 
matter of computational convenience. From physical considerations it is clear that the 
boundary should be chosen neither close to the crack nor close to the support point, be- 
cause of the stress-field disturbances near these positions. Preliminary studies estab- 
lished that V/W equal to 1. 5 was about optimum and that minor variations from 1. 5 had 
a negligible effect on the K calibration. A more detailed discussion of the effect of the 
choice of V/W has been given for the case of pure bending in reference 4. 

For each set of selected values of the primary variable a/W and the parameters 
S/W and 2y Q /W, the collocation computation was carried out four times, using succes- 
sively 15, 18, 21, and 24 boundary stations. In no case was the variation among the 
four d^ values so obtained as great as 1 percent, and in most cases it was a small 
fraction of 1 percent. From this, together with the nature of the trends of the dj values, 
it was concluded that those values corresponding to 24 boundary stations were very close 
to the limit values for large numbers of boundary stations. These values were accord- 
ingly used to calculate the results reported herein. 

The d-. values so obtained were used to calculate values of the square of the dimen- 

A 2 2 2 3 2 

sionless stress-intensity coefficient, namely Y = K B W /M , where M = PS/4. In 

this general form the results are of more immediate utility than they would be in the form 

of values of K/M for a specimen of specific, arbitrary dimensions. 


RESULTS AND DISCUSSION 

Dependence of on S/W and 2y 0 /W 

2 

Figure 2 shows plots of Y against S/W for constant a/W values of 0. 1, 0.2, and 
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0. 3. In each case the corresponding value of Y for pure bending is shown as a horizon- 
tal dashed line for comparison. The values for pure bending were taken from reference 4 
and, of course, are independent of S/W. It is clear that the value of Y for three-point 
loading is always lower than that for pure bending, but that the difference decreases with 
increasing S/W. For a/W equal to 0. 3, the ratio of the value of Y for three-point 
bending to that for pure bending is 0. 87 when S/W is equal to 4 and increases to 0. 94 
when S/W is equal to 10. 

2 

The trend of the curves in figure 2 indicates that Y becomes increasingly sensitive 
to S/W as this ratio decreases. This is a consequence of the increasing complexity of 
the overall stress-field pattern with decreasing S/W and is the reason why computa- 
tions were not conducted for S/W less than 4. It was considered that a K calibration 
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Figure 2. - Dependence of square of stress-intensity coeffi- 
cient for three-point bending on ratio of support span to 
specimen depth. 


parameter 2y Q /W on the K calibration can 


for S/W substantially less than 4 would 

be of very dubious accuracy, in fact, 

more misleading than useful. 

o 

The values of Y plotted in figure 2 

are those for the upper limit values of the 

parameter 2y Q /W, equal to 

2(1 - a/W) 2 W/3S. If the corresponding 

values for the lower limit of 2y Q /W, 

equal to zero, had also been plotted, they 

would have been virtually indistinguishable 

from their companions. For S/W equal 

2 

to 4, the values of Y for the lower 

limit of 2y Q /W were about 1 percent 

greater than the plotted values. For S/W 

o 

equal to 8, the values of Y for the lower 
limit were less than 1/2 percent greater 
than the plotted values. For practical 
purposes, therefore, the effect of the 
3 considered negligible. 


CALIBRATIONS OF K FOR S/W EQUAL TO 4 AND 8 

Figure 3 shows plots of the computed values of Y against a/W for S/W equal to 
8 and 4, again for the upper limit values of the parameter 2y Q /W. 

The following empirical equations are compact expressions of the same results for 
values of a/W up to 0. 35 

Yg = 134 a/W - 247 (a/W) 2 + 813(a/W) 3 (4) 


Y 2 = 130 a/W - 262 (a/W) 2 + 820(a/W) 3 


(5) 


2 o 

where Y g refers to S/W equal to 8, and Y^ to S/W equal to 4. These equations 
were obtained by least- squares-best- fit computer programs for cubics in a/W, incor- 
porating the known condition that K should be zero when a/W is zero. Only the results 
for a/W up to 0. 35 were used in fitting the equations since it is neither desirable nor 
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necessary to use bend specimens having cracks deeper than about 0. 35 W. The results 
are fitted by equations (4) and (5) with an average deviation of less than 0. 2 percent, the 
maximum deviation not exceeding 0. 5 percent. The corresponding equation for pure 
bending (ref. 4) is 


y| b = 139 a/W - 221(a/W) 2 + 783(a/W) 3 (6) 

Comparison With Results of Independent Studies 

The results of the present study for S/W equal to 8 are compared with the results 
of references 5, 6, and 7 in table I (p. 3). In general, the present results are in between 
the experimentally derived results of references 5 and 6, those of reference 5 being 
lower than the collocation values and those of reference 6 being higher. The results from 
reference 6 in table I are calculated from the fitting equation to the experimental data 
that is given in that reference. Values for a/W greater then 0. 35 are omitted because 
the fitting equation obviously does not represent the experimental results in this range. 

While the present results and those of reference 6 agree fairly well, the agreement 
is not as good as that demonstrated in reference 4 between the boundary collocation re- 
sults for pure bending and experimental results for four-point bending. Furthermore, 
the agreement between the two sets of experimental results for three-point bending 
(refs. 5 and 6) is distinctly poorer than the agreement of either with the present boundary 
collocation results. These facts would appear to indicate that it is inherently more dif- 
ficult to obtain an accurate K calibration for three-point bending than it is for four-point 
bending. 

Reference 7 cites only three actual results from Bueckner's unpublished analysis of 

o 

notched bars in bending. These results, converted to values of Y , are listed in table I. 
For a/W equal to 0.2, the agreement with the present result is quite good, in fact, 
better than the agreement with either set of experimentally derived results. The agree- 
ment for a/W values of 0. 1 and 0. 5, however, is no more than fair. This is a matter 
of academic rather than practical concern since the optimum range of a/W for Kj c 
testing with bend specimens is between 0. 15 and 0. 35. 


K Ic Measurement Capacity in Relation to W 

As discussed at some length in reference 1, the greatest value of Kj c that can be 
measured accurately with a given bend specimen, or the Kj c measurement capacity of 
the specimen, depends on both the depth W and the thickness B of the specimen. The 
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Relative crack depth, a/W 

Figure 4. - Dependence of plane strain crack 
toughness measurement capacity on relative 
crack depth for single-edge- notch specimen 
subjected to various types of loading. 


dependence on B is unrelated to the pres- 
ent work and will not be discussed here. 
Given a bend specimen of adequate thick- 
ness, the symbol Cj^- is used to denote 
the maximum value of Kj c that can be 
measured with acceptable accuracy by 
using that specimen. The currently ac- 
cepted criterion for evaluating C-^ is 
that the nominal stress at the position of 
the crack tip should not exceed the yield 


strength of the material in a valid Kj c 
test; that is, 6M/B(W - a) 2 should not 


exceed cryg. 


Substituting in the expression 


Y 2 = K 2 B 2 W 2 /M 2 , and transposing give 


C nC /<T YS W = y2(1 " a /W) 4 /36 (7) 

2 2 

The dimensionless quantity C^/ CygW is a ^Ic measurement efficiency factor which 
can be similarly evaluated for other types of specimens, as in reference 1. In all cases 
it is a function of a/W since Y is a function of a/W. In the special case of three- 
point bending it is also somewhat dependent on S/W, to the same extent that Y depends 
upon S/W. For any given values of a/W and S/W, there is a unique value of c ^/ a ygW, 
and therefore, is proportional to the yield strength of the material and to the square 

root of the specimen depth W. 

9 o 

Plots of Cj^/OygW against a/W for three-point bending with S/W equal to 4 and 8 
are shown in figure 4. Also shown for comparison are similar plots for pure bending and 
uniform tension, taken from references 4 and 3, respectively. The measurement efficien- 
cy is greatest when a/W is about 0. 25 in all cases. Furthermore, over the range of 
a/W from 0. 15 to 0. 35 the measurement efficiency is within 10 percent of the maximum 
value in each case. The differences in Kj c measurement capacity according to the dif- 
ferent methods of loading shown in figure 4 are insufficient to be of much practical impor- 
tance. The degree of uncertainty about the criterion used to calculate the measurement 
efficiency factors precludes drawing any fine distinctions in this respect. For practical 
purposes, therefore, it is reasonable to assume that the measurement capacity of a 
single -edge -notch specimen will be independent of the manner in which it is loaded. A 
convenient working rule for all single -edge -notch specimens is that Cj^- is about 
aYsW^ 2 /2 (for a/W in the range 0. 15 to 0. 35). 
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CONCLUDING REMARKS 


In this report the results obtained by the boundary collocation procedure are expressed 
in general form in terms of the square of the dimensionless stress-intensity coefficient 
Y 2 = K 2 B 2 W 3 /M 2 . In the case of pure bending, is a function of the relative crack 

p 

length a/W only. In the case of three-point bending, the computed values of Y were 
appreciably lower than the corresponding values for pure bending (computed previously by 
the same method), the more so the smaller the ratio of support-span to specimen depth 
S/W. Thus, a different K calibration plot of Y against a/W is needed for each dif- 
ferent value of S/W that is used in plane-strain crack toughness testing with three-point 
bend specimens. Accurate plots of Y against a/W and fitting equations are given for 
S/W equal to 8 and 4. 

p 

For S/W equal to 8 the maximum deviation of Y for three-point loading from that 

for pure bending was about 7 percent. For S/W equal to 4 the maximum deviation was 

2 2 

about 13 percent. Since Y is proportional to K , the corresponding deviations for K 

are approximately half as great. It is considered that results for S/W substantially less 

than 4 would be of dubious accuracy and probably more misleading than useful. 

The K calibration relation of Y to a/W is also slightly affected by the spread of 

the contact pressure region around the center loading point. This spread depends upon 

t he yield strength and toughness of the material and the size of the specimen tested. Over 

the practical range for acceptable plane-strain crack toughness measurements, the effect 

of this factor was at most 1 percent and can therefore be considered negligible. 

The results agree fairly well with published results derived from experimental com- 

o 

pliance measurements; one set gave higher values of Y than the present method, and 
the other set gave lower values. 

The plane-strain crack toughness measurement capacity as related to the depth of 
three-point bend specimens was estimated to be somewhat lower than that of four-point 
bend specimens, but for practical purposes this difference is probably of little importance 
For all single-edge-notch specimens, whether tested in tension or bending, the measure- 
ment capacity is greatest in the range of a/W between 0. 15 and 0. 35. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 20, 1965. 
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APPENDIX A 


SYMBOLS 


a 

crack depth 

r 

B 

specimen thickness 


C IK 

estimated maximum value 

S 

of K t that can be mea- 
Ic 

sured with specimen of 
given dimensions and 

V 


yield strength 

W 

d 2n, d 2n-l 

coefficients of Williams 
stress functions 

x,y 

K 

stress-intensity factor of 
elastic-stress field in 
vicinity of crack tip 

Y 

K lc 

critical value of K, at 
point of instability of 
crack extension in first 
or open mode, a meas- 
ure of plane-strain crack 

2y 

J o 


toughness of material 

e 

M 

applied bending moment 

<7 

m 

number of selected bound- 
ary stations used in col- 

X 


location computation 

a 

y 

P 

total load applied to spec- 



imen 

ct ys 

T xy 


x 


polar coordinate referred 
to crack tip 


span 

distance from crack to 
boundary selected for 
collocation analysis 

specimen depth 

Cartesian coordinates 
referred to crack tip 


dimensionless stress- 
intensity coefficient, 

kbw 3/ vm 


length over which applied 
load at center was 
assumed to be distri- 
buted 


polar coordinate referred 
to crack tip 

stress component in x- 
direction 


stress component in in- 
direction 

0. 2 percent offset tensile 
yield strength 

shearing stress component 

stress function 
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